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Abstract 

Using tadpole inproved actions we investigate the consistency between different 
methods of measuring the aspect ratio renormalization of anisotropic-lattice gluons 
for bare aspect ratios xo = 4, 6, 10 and inverse lattice spacing in the range a^ 1 — 
660 — 840 MeV. The tadpole corrections to the action, which are established self- 
consistently, are defined for two cases, mean link tadpoles in Landau gauge and 
gauge invariant mean plaquette tadpoles. Parameters in the latter case exhibited 
no dependence on the spatial lattice size, L, while in the former, parameters showed 
only a weak dependence on L easily extrapolated to L = oo . 

The renormalized anisotropy xb, was measured using both the torelon dispersion 
relation and the sideways potential method. We found good agreement between 
these different approaches. Any discrepancy was at worst 3—4% which is consistent 
with the effect of lattice artifacts that for the torelon we estimate as 0{asa 2 s / R 2 ) 
where R is the flux-tube radius. 

We also present some new data that suggests that rotational invariance is es- 
tablished more accurately for the mean-link action than the plaquette action. 
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1 Introduction 



In principle, an improved action makes it possible to achieve lattice volumes large 
enough to overcome finite size effects at a computational effort low enough to obtain 
measurements with good statistical errors. However typical masses for heavy states 
may be similar to or larger than the inverse lattice spacing on the associated coarse 
lattice. Propagators for such states then decay too fast in lattice units for accurate 
measurement. 

The problem can be overcome by tuning couplings associated with the "time" di- 
rection so that the temporal lattice spacing a% is much smaller than the spatial lattice 
spacing a s . Refined measurements are then possible while retaining the computational 
advantages of the improved action and coarse spatial lattice. Such anisotropic actions 
have already been successfully applied to the glueball spectrum p], |2|, the spectrum 
of excitations of the inter-quark potential ||, heavy hybrids j|, ||, and the fine struc- 
ture of the quarkonium spectrum Q. They are also expected to be a powerful tool in 
extracting excited-state signals, obtaining high-momentum form factors, pushing ther- 
modynamic calculations to higher temperatures, and in the calculation of transport 
coefficients. 

In this paper, we consider one formulation of the anisotropic action for the pure 
Yang- Mills sector |], |] : 



-l J 5 P s , s i 1 R ss ,s' 1 Rs's',s\ osr^ /4 Ps,t 1 Rss,t 



9 \3 uj 12 ui 12 u{ J P h X °\3^ Ylu\u\ 



x,s>s' v 8 s s ' x,s 

Here s, s' run over spatial directions, P s>s ' is a 1 x 1 plaquette and R s > s ',s is a 2 x 
1 rectangle. The coefficients of these terms are chosen so the action has no 0{a 2 s ) 
discretization errors in tree-level pertubation theory, xo is the anisotropy parameter 
and is equal to the aspect ratio of the spatial and temporal lattice spacings, a s and at 
at tree-level. At higher orders in the perturbative expansion, this aspect ratio recieves 
quantum corrections, so a renormalized anisotropy determined from a physical process, 
Xr, differs from xo a t 0(a s ). Tadpole improvement (TI) of the perturbative expansion 
is achieved by tuning the input spatial and temporal link parameters u s and ut for 
self-consistency at each choice of (3 and xo- 

Here we report on the accurate determination of the tadpole parameters in the 
plaquette and Landau mean-link formulations. (Non-tadpole-improved actions have 
been studied in ref. || [Toft.) We compare methods for measuring the renormalized 
anisotropy and discuss the importance of lattice artifacts. In order that the anisotropic 
formulation can be used with confidence, physically distinct methods for determining 
the renormalized anisotropy, \R = a s/ a ti should agree. Discrepancies in the results for 
Xr give a measure of lattice artifacts and an indication of the importance of further 
improvements. Our results have been used in applications of anisotropic lattices to 
heavy quark hybrid states ||. 

In section ^ we report on the determination of the tadpole parameters; in section 
|3| the renormalized anisotropy measurements are presented, in section || we discuss the 
results and present conclusions. 
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2 Tadpole Parameters 



In both the plaquette and Landau mean-link formulations, the tadpole parameters are 
determined self-consistently and are defined by 



(P,. 



>V4 



s I Landau j 



Ut 



1 



t Landau 



plaquette, 



Landau, 



(2.1) 



where the Landau gauge is defined by the field configuration which maximizes [11] 



u,, af, 

x a M a 



Un(x)UJx + p,) 



(2.2) 



with respect to gauge transformations. We denote the gauge coupling in the two 
schemes as (5l and (3p respectively. 

2.1 Plaquette Tadpoles 

In the plaquette scheme of Eqn. ( |2.lD , ut = 1 and u s is determined self-consistently. The 
expectation value of the spatial plaquette is computed for a range of input parameters 
u s close to and spanning the self-consistent value, u*. A linear interpolation is sufficient 
to give an accurate value of u*. This value is then checked in a further Monte-Carlo 
simulation. The values of ti* for a set of couplings (3p and tree-level anisotropics xo 
are given in Table |l[ For the lattice sizes used in simulation, the plaquette expectation 
value is found to be independent of the volume at the four-significant-figure level. An 
extrapolation to infinite volume is unnecessary. 





Xo 


K 




(P SyS >) at u* s 


2.1 


6 


0.774166 


0.3592 


0.359136(72) 


2.3 


4 


0.793925 


0.3973 


0.397267(52) 


2.3 


6 


0.791062 


0.3916 


0.391621(73) 


2.3 


10 


0.789695 


0.3889 


0.388650(85) 



Table 1: Self-consistent plaquette tadpole parameters. 



2.2 Landau Gauge Fixing 

The maximization of F(U2) in Eqn. ( |2.2j ), where US = g{x)U a {x)g\x + /j,) with re- 
spect to a gauge transformation {g(x)}, was carried out using the conjugate-gradient 
method modified to deal with the group structure of the link elements. At each stage 
of the iteration the the appropriate conjugate-gradient vector {v(x)}, is computed as 
a covariant derivative [12]. 



d_ 

drj(x) 



< x ) = T^ F ( U P\v=o > (2-2-1) 



where g(x) = e n ^' T and T are the generators of SU(3) . The vector {v(x)} has 
elements lying in the Lie algebra of SU(3) . For each vector a series of group elements 
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is constructed in the associated one-parameter subgroup, each element being a given 
power of the preceding one: 

5i = exp(e^) , g n = ( 5n _i) p , 1< n < N . (2.2.2) 

The local maximum in the direction v is found by evaluating F(Uji p ) in descending 
order from p = N to p = 1 . The value of e, at any stage, is reduced until either the 
sequence exhibits a maximum or remains constant within a preset tolerance. 

For large xo t ne terms including Ut dominate the expression for F and near to 
the maximum relatively small changes in F correspond to quite large changes in u s . 
Consequently, to be sure that u s as well as ut is accurately calculated, the maximum 
must be found to a sufficiently high accuracy. The criterion chosen was that 5F/F < 
2 x 10~ 6 where 5F was the accumulated change in F for three consecutive iterations 
of the conjugate-gradient algorithm for which the individual changes in F were non- 
zero. This criterion accounts for the observation that from time to time close to the 
maximum the change in F was zero. We chose e = 2 x 10 -5 , p = 4, and N = 6. 



2.3 Landau Gauge Fixed Tadpoles 

The self-consistent value of u = (u s ,ut) can be determined by a generalized Newton- 
Rapheson method or by linear interpolation. The latter method was used for the mean- 
link tadpoles and was implemented by choosing four input tadpoles Ui, i = 1, . . . ,4 
and measuring the four corresponding output tadpoles u^, i = 1,4. A linear map is 
assumed for the incremental vectors r$ = (itj — u±), — ► r[ = (u^ — u[), i = 2, 3, 4: 

r» = MrJ , i = 2,3,4 , (2.3.1) 

where M is a 2 x 2 matrix that is determined from the the measured images of r2 , r3 
and checked for consistency against the measured image of r3 . The self-consistent 
tadpole u* is then predicted to be 

u* = Ul + (l-M)-Vi-ri) ■ (2.3.2) 

The self-consistency of u* is then checked computationally. This method was found to 
be very reliable and only subject to minor adjustments to account for statistical errors. 

Moreover, a meaningful statistical error can be assigned to u* from the image under 
M _1 of the statistical error box deduced for u* from the measurements values of the 
tadpole parameters. Let the image of u* under the map M be u*' with statistical 
error 5u*'. The value u* is acceptable if in terms of components u*' — 5u*' < u* < 
u* f + 5u*', and the error quoted on u* is 5u* = M _1 Su*'. A typical example is for 
/3l = 1-8, xo = 4 on a 6 3 x 24 lattice where 

/ 0.341 -0.227 \ . *, / 3xl0~ 4 \ . , / 1 x 10~ 4 \ 
Ml = { _ .033 0.957 J' 5U = ^3x10-5 J =* 6U = (2x10-5 J" 

(2.3.3) 

In practice it is found that the error on u* is typically reduced by a factor of 3 compared 
with the statistical error from the simulation verifying the self-consistency. 

The simulation for the mean-link tadpoles was done on the Hitachi SR2201 comput- 
ers at the Cambridge High Performance Computing Facility and the Tokyo University 
Computer Centre. The lattices used were L 3 x T with T = xoL and L = 6, 8 in all cases 
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except one. The high accuracy demanded by the maximization process was very time 
consuming. Also, many configurations were required to achieve the desired statistical 
accuracy. Consequently only one example with L = 10 was done as a check on the 
finite-size scaling ansatz, u*(L) = u*(oo) + A/L 2 . Typically for 10 3 x 40 about 400 
conjugate-gradient iterations were needed taking about 60 seconds per iteration. In 
the cases L = 6, 8, 10 for /3l = 1-8 xo = 4, the fit to the finite-size scaling ansatz was 
very good as can be seen in Figs, [l] and ||. For other cases the linear extrapolation was 
assumed to hold. The values measured for u*(L) and the extrapolation to L = oo are 
given in table ||. 





Xo 


L 


4 


6 


8 


10 


oo 


1.8 


4 


0.7165(2) 
0.98124(3) 


0.7244(1) 
0.98201(2) 


0.7260(1) 
0.98222(2) 


0.7266(1) 
0.98227(2) 


0.7279(2) 
0.98243(3) 


1.8 


6 


0.7115(2) 
0.99167(3) 


0.7191(1) 
0.99194(2) 


0.7202(1) 
0.99200(1) 




0.7216(3) 
0.99208(3) 


1.7 


4 




0.7127(1) 
0.98105(2) 


0.7143(1) 
0.98188(2) 




0.7164(3) 
0.98295(5) 


1.7 


6 




0.7075(1) 
0.99149(1) 


0.7086(1) 
0.99153(1) 




0.7100(3) 
0.99158(3) 



Table 2: Self-consistent Landau mean-link tadpole parameters for various lattice sizes, L, 
including the 1/L 2 extrapolation to L — oo. 

3 Anisotropy 

In this section we report on two methods for determining \R- The first uses the disper- 
sion relation for the torelon (jDjI and the second uses the comparison of the potential 
measured in the fine and coarse direction using Wilson loops. 

3.1 The Torelon 

The lattice considered is S 2 x LxT with typical values S = 8 in the x and y directions, 
L = 3, 4, 5 in the z direction, and T = 50 in the fine or t direction. The lattice has 
periodic boundary conditions and the torelon is created by a Polyakov line that loops 
around the lattice in the z direction. The Polyakov line is associated with a particular 
point in the (x, y)-plane and is constructed from links which are covariantly smeared 
in a manner similar to APE smearing: 

z=L-l 

T(x,y,t) = Tr J] W z (x,y,z,t) , 

z=0 

l 2 D 2 

W z (x,y,z,t) = (l + —) m U z (x,y,z,t), (3.1.1) 

where D is the appropriate covariant derivative. Typically, I = 1.5, m = 10. The state 
with momentum p = (p x ,p y ) = (n x ,n y )(2ir / Sa s ) is 

T(p, t)=J2 T (z> 2/> t) e iip * x+p y y) . (3.1.2) 
x,y 
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The torelon propagator 



-. t'=T-l 

G T (p,t) = - (T(p,t')T*(p,t + t')) , (3.1.3) 

t'=0 

is measured for various choices of momentum p and a correlated simultaneous fit using 
SVD decomposition is made to the relativistic dispersion formula: 

G T (p,t) = c(p) e- E ^ , 
aJp 2 + M 2 

E{p)a t = — V 1, (3.1.4) 

XR 

where My is the torelon mass and i = i/at- The momenta used were n 2 = n 2 + 
= 0,1,2,4,5 with 5 = 8. An example of the fit obtained is shown in Fig. |3] for 
Pl = 1-8, Xo = 6. The fit is over the range 2 < t < 12 and has xV-^df = 0.96 for 
43 d.o.f. The fit is good and gives (M T a s ) 2 = 1.85(4), X R = 3.61(2). In Fig. | the 
dispersion curves E(p 2 ) versus n 2 are plotted for 0l = 1-8, xo = 4, L = 3,4,5. From 
both Figs. H and [| it is clear that rotational invariance is established in the coarse xy 
directions. The good simultaneous fit to the torelon propagators in all cases shows that 
the continuum dispersion is well satisfied for momenta at any angle to the coordinate 
axes. 







L 




Xo 


3 


4 


5 






XR 


{M T a s y 


XR 


{M T a s f 


XR 


(M T a s y z 


1.8 


4 


3.61(1) 


0.683(1) 


3.61(3) 


1.85(4) 


3.57(5) 


3.39(13) 


1.8 


6 


5.32(2) 


0.502(1) 


5.28(3) 


1.49(3) 


5.28(5) 


2.84(9) 


1.7 


4 


3.56(2) 


1.16(2) 


3.56(4) 


2.86(9) 


3.66(9) 


5.51(36) 


1.7 


6 


5.28(2) 


0.89(2) 


5.28(5) 


2.36(6) 


5.26(6) 


4.25(13) 



Table 3: Renormalized anisotropics and dimensionless torelon masses for various values of 0l 
and xo from a simultaneous fit on 8 2 x L x 50 lattice. 

In table ^ we give the full set of results for the torelon anisotropics we have measured. 
There is no noticeable dependence within errors of XR on L an d the dependence on 
is small for the two values of /3 used. 

From the values Mr(La s ) we can determine the string tension, a in units of a~ 2 . 
We assume that Mt(Lo, s ) = a(La s )La s where a(La s ) is the string tension modified by 
finite-size corrections |l4], [l5|]: a{La s ) = a + D / (La s ) 2 . A plot of a{La s )a 2 versus 1/L 2 
is shown in Fig. |B| with a linear fit which has xV-^df = 0.05 with aa 2 = 0.394(3). In 
another calculation we have determined the absolute value of at in MeV by computing 
the splitting AMps = Mx(lP) — Mx(lfi') for the T meson system using 0(mv e ) 
NRQCD |i"6| . In table | we give a^ 1 , a' 1 , aa 2 sl the coefficient D, and / l\Mps- 



This latter ratio is experimentally close to unity but in ref. 16] using isotropic /3 = 6.2 
UKQCD configurations this ratio was found to be about 1.25 which, from table ||, 
agrees with our finding except possibly for the (/3 = 1.7, xo) lattice which has the 
most coarse lattice spacings. It is generally accepted that the discrepancy is due to 
quenching and the significant outcome is that our results from the anisotropic lattice 
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agree well with those of a spatially finer isotropic lattice of a = 0.06fm, indicating that 
we have correctly reproduced the physics expected for this comparison on lattices with 
a s = 0.25 — 0.30fm. A naive estimate for the coefficient D is obtained by calculating 
the contribution to M{L) from the zero-point fluctuations of the torelon regarded as 
a periodic flux tube. The outcome is D = — it/3 ~ —1.05 which is compatible with 
our fit of D ~ —1.35(5). The same calculation predicts that there is no L-independent 
constant contribution to M(L) and this is verified by our fits. 





Xo 


M b 


a f -1 MeV 


aJ x MeV 


aa 2 s 




D 


1.8 


4 


5.32 


2876(75) 


797(21) 


0.422(7) 


~1.18(3) 


-1.31(6) 


1.8 


6 


4.75 


4503(45) 


839(9) 


0.394(3) 


1.19(2) 


-1.42(3) 


1.7 


4 


5.97 


2353(38) 


661(11) 


0.513(12) 


1.08(2) 


-1.39(13) 


1.7 


6 


5.56 


4112(150) 


779(28) 


0.469(8) 


1.21(5) 


-1.38(8) 



Table 4: The inverse lattice spacings for mean- link tadpole improved actions from the 0(mv 6 ) 
NRQCD measurement of AMps = Mr(lP) — Mr(lS') to determine a^ 1 and using the torelon 
anisotropy to infer aj 1 . The string tension is found from the L— dependence of the torelon 
mass using finite-size scaling ansatz a{La s ) = a + D / (La s ) 2 . 



The torelon dispersion relation was also studied on the plaquette-tuned parameters. 
For these simulations, lattices of extent (8 x 10) x 4 x Nf were used, to enable us to 
investigate a wider range of momenta combinations. At momenta close to the cut-off, 
we expect the discretization errors to be 0(aj.p 4: ,a s d 2 .p 2 ), thus in order not to con- 
taminate our determination of \R with large discretization errors, we first determined 
the range of momenta over which a good correlated fit to the continuum dispersion 
relation could be made. The data for different momentum ranges from simulations at 
f3p = 2.1, xo = 6 were tested and these results are shown in Table ||. Both the measured 
anisotropy renormalizations and the torelon rest energies determined are all consistent 
within statistical precision up to the highest momentum measured. For the subsequent 
computations of the anisotropy, presented in Table ^, the largest momentum used was 
the (1,1) data (shown in bold in Table |5|). 



Maximum 


xr/xo 


a t Mx 


xV^df 


(Px,Py) 








(1,1) 


0.949(25) 


0.2884(17) 


1.19 


(0,2) 


0.972(23) 


0.2886(16) 


0.99 


(1,2) 


0.968(17) 


0.2885(15) 


0.85 


(2,1) 


0.976(14) 


0.2885(15) 


0.89 


(2,2) 


0.970(11) 


0.2883(14) 


0.96 



Table 5: Dependence of the renormalized anisotropics and torelon masses on the 
highest momentum used in the fit to Eqn. ( |3.1.4 ) at j3p = 2.1, xo = 6. 



3.2 The Sideways Potential 

In the sideways potential method |9j, O], a coarse direction, z, on the anisotropic lat- 
tice is chosen to be the time direction. There are then two types of spacelike direction, 
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Pp 


Xo 


Xr/Xo 


X 2 /N d{ 


2.3 


4 


0.941(13) 


1.04 


2.3 


6 


0.954(23) 


0.74 


2.3 


10 


0.940(13) 


1.41 


2.1 


6 


0.949(25) 


1.19 



Table 6: Renormalized anisotropies for various values of f3p and xo using plaquette-tuned 
mean-link parameters. 



coarse and fine. The measurement of the potential between static quarks with a sep- 
aration lying in the plane of coarse links is compared with the measurement of the 
potential when the separation lies along the line of fine links. The demand that the two 
measurements yield the same function of physical distance determines the renormalized 
anisotropy. Points in the coarse-coarse plane are denoted by x = (x, y) and points in 
the fine direction by t where x, y, z, t are integers. 

We measure appropriate spatial Wilson loops W ss (x,z) and also loops using the 
fine direction Wt s (t,z). We define 

T/ ,- v , ( W ss (x,z) \ ( W ts (t,z) \ 

V 8 (x, z) = log — , V t (t, z) = log — — — — - . (3.2.1 

\W ss (x,z + l)J \W ts {t,z + 1)J 

As z — > oo, V s (x,z) — > V s (|x|), and V t (t,z) — > V t (t) where V^(|x|) and V t (t) are the two 
versions of the interquark potential. For a physical distance r we have \x\a s = tat = r ■ 
We therefore estimate the renormalized anisotropy xr by tuning it so that = 
Vt{t/xR)i where the right side is evaluated by means of linear interpolation between 
the values measured at integral t. It is implicit in the method that there is effective 
rotational invariance in the x-plane. We find that if we exclude potentials at the smallest 
distances, |x| = 1, y/2, then the values of xr associated with different directions in the 
x-plane generally agree within errors. The agreement is particularly good if the links 
are smeared in an appropriate manner. The results for the renormalized anisotropy, 
Xr, for both mean-link and plaquette schemes are shown in table f?|. 

An alternative approach to making the comparison is to fit the measured potentials 
with the forms 

a s V s (x) = a s V + cat x + ~ 

x 

a s Vt(t) = a s V + aa s a t t + . (3.2.2) 

a t t 

The renormalized anisotropy, XRi 1S then determined from the ratio of the coefficients 
of the linear terms in the two cases, xr can m principle be determined from the ratio 
of the coefficients of the coulombic terms however such an estimate depends on short 
distance effects and is inherently more sensitive to discretization errors. 

This approach was tested on an 8 3 x 48 lattice at j3p = 2.3, %o = 6. Prior to 
measurement, the lattice is blocked to an 8 volume by thinning the time-slices of spatial 
links and replacing the temporal links with the product of six fine links connecting the 
appropriate time-slices. The t and z axes are then interchanged, since as before we 
want to use a coarse direction for Euclidean decay. Then a standard APE-smearing 
algorithm is applied to links on the new time slices. These degrees of freedom are then 
used to measure the potential in the (x, y) and z axes in the standard fashion. Fig. || 
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mean-link 


plaquette 


P 


Xo 


XR 


P 


Xo 


XR 


1.8 


4 


3.67(2) 


2.1 


6 


5.80(5) 


1.8 


6 


5.48(3) 


2.3 


4 


3.95(3) 


1.7 


4 


3.67(3) 


2.3 


6 


5.93(5) 


1.7 


6 


5.40(5) 


2.3 


10 


10.00(20) 



Table 7: Results for the renormalizcd anisotropy, xr, using the sideways potential method for 
both mean-link and plaquette tadpole-improved schemes 



shows the results of this simulation. The best fits of this data to Eqn. ( 3.2.2 ) for a range 
of static-source separations were computed and the string tension results are shown in 
Table ||. In Fit C, the coulomb term was poorly resolved and was thus fixed to zero. 
The use of these different ranges allowed us to investigate the systematic uncertainty 
and discretization errors in determining xr from this method. The three results are 
consistent within 2%. The anisotropy measurement from the full range of separations 
(Fit A, shown in bold in Table ||) is in good agreement with the torelon dispersion 
result of Table H. 



Fit 


Separations 


aa 2 s 


xV^df 


6aa s at 


xV^df 


xr/xo 


A 


1-7 


0.3293(33) 


1.13 


0.3492(26) 


1.21 


0.9430(12) 


B 


1-4 


0.3315(60) 


0.85 


0.3526(64) 


1.14 


0.9402(24) 


C 


3-7 


0.3506(32) 


0.70 


0.3664(34) 


1.11 


0.9569(12) 



Table 8: Results for the renormalized anisotropy, XR from fits to the static potential. 
In fit C, the Coulomb term was set to 0. 



Both the above approaches yielded results reasonably consistent with each other 
and with the toleron results at the 3% level. Discrepancies can easily be explained 
by the presence of discretization errors. Because the flux tube generating V^(|x|) has 
a coarse-fine cross-section and that generating Vt(t) a coarse-coarse cross-section we 
anticipate, as in the case of the toleron (eq.([4.5|), that some discretization error in XR 
that is 0((a s a s /R) 2 ) will remain. Since a s ~ 0.3 fm, R ~ 0.7 to 1 fm, and a s ~ 0.3, 
we expect discretization errors to be around 5%. 

Because we are concerned to make the long distance behavior consistent in both the 
fine and coarse directions it is advantageous to use Wilson loops of the largest possible 
spatial extent. However, in practice, the statistical errors on large Wilson loops grow 
exponentially with separation. In order to achieve acceptable errors we were restricted 
to using loops of size 2 or 3 in coarse lattice units. This is to be compared with the 
torelon, where the practical size is 3 or 4 lattice units. 

4 Discretization Errors 

One way of viewing discretization errors is to regard them as arising from the ab- 
sence of the correction term in the lattice action that yields continuum results. The 
corresponding term in the action density can be presumed to be a local (redundant) 
operator. The locality of the operator implies that while there may be a correction to 



9 



the mass per unit length of the torelon, there is no correction that directly alters the 
asymptotic proportionality of the torelon mass and its length. The appropriate length 
scale against which to measure lattice artifacts therefore is the radius R of the toleron 
cross-section or equivalently, Aq^ d or yfa where a is the string tension. We expect 
therefore that the 0{a 2 s ) errors in \R are proportional to (a s /R) 2 . 

A persuasive plausibility argument for this behavior is as follows. The torelon 
state of momentum p and mass Mr is an eigenstate of the transfer matrix T with an 
eigenfunctional labelled by p and Mr and eigenvalue exp(— E(p 2 )at/xR)- Suppose an 
operator is added to the action which corrects for the 0(a s a 2 ) errors. This operator 
will be local on the scale of the lattice and the effect on E([p 2 ) can be estimated using 
first-order perturbation theory: 

6 a_Mp 2 )_ = asf{p2jMTh (41) 
XR 

where p 2 = p 2 a 2 , Mt = Mra s , an d /(p 2 ,Mr) is the dimensionless matrix element of 
the added operator which is proportional to a 2 by construction. The question is what 
scale balances the dimension of a 2 ? We find 



"s ■ 



a s <f>'(0,M T ) , (4.2) 



Xr 

= a s (<ft(Q,M T )-M T (ft'(Q,M T )) , (4.3) 

%XR 

where (ft(p 2 , M) = a s E(p 2 ) f (p 2 , Mt) and (ft 1 = d(ft/d(p 2 ). It must be that the change in 
the action by a local operator corresponds to a change in the string tension, a. Hence, 
since My = aL, it follows that 5M T oc M T as this is the only parameter depending on 
L, and from eqn. ( |4.3| ) this implies that 

0(0, M T ) ~ Mf > 0'(O,M T ) ~ Mr-independent constant . (4.4) 



It is conceivable that an accidental cancellation between the two terms in eqn. (|4.3[) 
would allow different behavior to be inferred for (ft and (ft 1 but these results must hold 
true for all possible local perturbations, not just the particular one that eliminates 
0(a s a 2 ) errors. We consider this kind of cancellation to be unlikely. Substituting this 
behavior for (ft 1 into Eqn. (|1.2j) we find that 

5xr r* a s (A <\ 

— R ~C T2 , (4.5) 

where R is a typical length scale associated with the torelon which cannot be M^ 1 , for 
example, the flux tube radius. Alternatively, R~ l can be taken to be y/a, A.qcd- 



5 Conclusions 

In this paper we have investigated various methods for measuring the renormalized 
aspect ratio for pure QCD on an anisotropic lattice. In the main we used a tadpole- 
improved action with the Landau gauge mean-field definition for the tadpole param- 
eters, but we also included results for the action with tadpoles defined by the mean 
plaquette. The object of the investigation was to assess the consistency of different 
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methods in order to judge the effectiveness of the improvement scheme. In principle, 
measurements of the anisotropy from different physical probes should agree close to the 
Euclidean-symmetric continuum limit. The bare anisotropy xo is renormalized by the 
effect of operators which are irrelevant in the neighborhood of the fixed point control- 
ling the continuum limit. From a renormalization group (RG) view point the location 
of this fixed point is ambiguous up to redefinitions of the RG transformation used to 
locate it. This ambiguity is due to redundant operators |l8| , |i~9| , 20 1 which have no effect 
on physical observables of the continuum theory. Consequently, continuum actions with 
different \R can differ only by redundant operators since they must correspond to the 
same continuum physics. Thus, while xr can be changed by tuning xo> an action with 
no lattice artifacts must give rotationally invariant physical results and, consequently, 
any physical method for measuring XRi such as the ratio of two physical observables, 
must give the same answer. In as much as this is not the case the differences will give 
an estimate for the effect of lattice artifacts and the necessity for improvement. 

The spatial lattice spacing a s was measured in a separate NRQCD simulation by 
fitting the IP — 15 mass-splitting AMps for bottomonium. This estimate can be 
compared with the string tension, deduced from a D/L 2 extrapolation to oo of the 
torelon mass. The coefficient D was found to be ~ —1.38 which is in tolerable agreement 
with —it/3 predicted from analysis of flux-tube fluctuations J15|| . The ratio ^/a/AMps 
is ~ 1.2 in agreement with earlier NRQCD analyses [21, |lq] . The departure from 
the experimentally observed value of ~ 1 is attributed to quenching. The values for 
both the mean-link and plaquette tadpoles were found self-consistently. This was very 
resource intensive in the mean-field case since it required a very accurate gauge fixing 
to Landau gauge and the self-consistent iteration is in the 2D space of (u s ,ut) which 
requires additional effort. Also, for given xo the tadpoles showed an L dependence 
and hence required an extrapolation to L = oo. This in turn requires accurate data 
for a good fit. In contrast, the plaquette tadpole was easily found and there was no 
discernible L-dependence. 

The torelon method deduces xr from a fully correlated fit to the dispersion relation. 
Statistical errors are produced by the fit. The sideways potential approach requires 
some method for matching the coarse and fine potentials and the error analysis is more 
complex because of both systematic and statistical errors and because the signal rapidly 
decreases as the loop size increases. Two methods were used to extract xr- The first 
compared the loop predictions for the coarse and fine potentials ( [3.2. 1 ) and deduced 
Xr from the rescaling of t needed for them to agree. This method was applied to 
actions with mean-field improvement. Typical loops were of edge length 2-3 in units of 
the coarse spacing a s . There is a clear difference in results compared with the torelon 
computation of about 3 — 4%. In the second approach the fine direction is first blocked 
by a scaling of xo to give a lattice that is approximately isotropic and then the potential 
is fitted to a standard form ( |3.2.2| ) and xr deduced by requiring that the linear slopes 
agree. The latter method has the advantage that it depends only on the long range 
structure of the potential and excludes the short range coulomb part. However, it does 
require large computing resources to extract a reliable signal at large separations, in this 
case up to 7 lattice units. This technique was applied mainly to actions with plaquette 
improvement. Although less extensively investigated, the results are consistent with 
the torelon computation. 

The observed discrepancies can be easily attributed to residual lattice artifact ef- 
fects. For the torelon we expect 5xr/xr ~ Ca s a 2 /R 2 which can be sizeable although 
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action 


V(lll) - VetiVS) 
V(2)-V(l) 


V(221) - V(3) 
U(3)-U(l) 


Wilson 


0.267(2) 


0.099(3) 


plaquette TI 


0.091(10) 


0.027(2) 


mean-link TI 


0.050(5) 


0.003(7) 



Table 9: Comparisons of the normalized deviation from rotational invariance between the 
Wilson action and the plaquette and mean-link tadpole-improved (TI) actions. The ratios 
shown are zero in the case where rotational invariance holds, and on these measures the mean- 
link TI action is superior. 



C is unknown. However, it should also be remarked that neither the torelon nor the 
second sideways potential method include the coulomb part of the static potential. 
Because of its short-range nature this part will be the most sensitive to the effect of 
lattice artifacts. In this case agreement of prediction for xr ls good. This may merely 
indicate that xr 1S being deduced from the properties of the flux tube in both cases and 
so they should agree in that artifact effects are the same. However, it is encouraging 
that different methods based on the long-range properties of the action are consistent 
and it is worth noting that the torelon flux tube was 3-5 lattice units in length whereas 
it was necessary to extend the potential method to 7 lattice units. Simulation times are 
correspondingly reduced for the torelon. Unfortunately, it was not possible to apply 
the second method to data gathered using the mean-field action because of limitations 
in resources. 

We do not find any evidence that the mean-link tuned action is superior to the 
plaquette tuned action in suppressing lattice artifacts, and it should be noted that the 
computation required to determine the mean-link tadpoles for the former case is very 
time consuming. However, there is other evidence that the mean-link action is superior. 
It is known to give smaller scaling errors in the NRQCD charmonium hyperfine splitting 
p2| , |23| ], better hadron mass scaling and to give a clover coefficient for SW fermions 
(with the Wilson gauge action) that agrees more closely [11] with the non-perturbatively 
determined value ]24] , pq ]. Using our configurations we present evidence that it also 
gives better rotational invariance for the static quark potential. In Fig. [7| we show 
the deviation of the potential from a fit to the standard form V(r) = b + c/r + or. 
Comparing the mean-link TI data with the plaquette TI data, it is clear that mean- 
link TI gives smaller deviations from the rotationally invariant fit. In table |9] we present 
some of the features of Fig. in a quantitative form. As expected, mean-link TI gives 
better rotational invariance at r = y/3, and it causes the two r = 3 potentials ((2, 2, 1) 
vs (3,0,0)) to agree better. 

As mentioned in the introduction, there are many applications for anisotropic lat- 
tices. It would be valuable to measure anisotropies for a wider variety of lattice spacings 
and bare anisotropies, using the methods we have investigated. It would also be very 
useful to calculate perturbative formulae for the anisotropy and lattice spacing as a 
function of the bare parameters of the action, (3 and xo H> 0> [2^] . Finally, it will be 
necessary to repeat this tuning process for improved light quark actions. Some of the 
requisite perturbative calculations have already been performed [27]. 
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Figure 1: Fit of u s (L) = u s (oo) + A s /L 2 for L = 6,8, 10. 
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Extrapolation of tadpole parameters 
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Figure 2: Fit of u t (L) = u s (oo) + A t /L 2 for L = 6, 8, 10. 
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Torelon propagators for various momenta 

P=1.8, % =4 




1 2 3 4 5 6 7 8 9 10 11 12 13 
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Figure 3: Fit of the relativistic dispersion Eqn. ( |3.1.4 ) to the torelon propagators for p 2 = 
n 2 (ir/A) 2 for n 2 = 0, 1, 2, 4, 5 corresponding respectively to the curves from the top downwards. 
The fit has x 2 /^Vdf = 0.96 with (M T a s ) 2 = 1.85(4), X R = 3.61(2). 
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Torelon dispersion on 8 xLx50 

M-8, % =6 
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1 2 3 4 5 

p 2 in units of (2rc/8) 2 



Figure 4: Plots of E(p 2 ) versus p 2 for L = 3, 4, 5 on 8 2 x L x 50 lattice. The lines show 
fit to E(p 2 ) =A + Bp 2 . 
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Figure 6: The potential at (3p = 2.3, \o = 6 between static sources separated along both 
spatial and temporal axes. 
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Figure 7: Deviation of static potential from fit, a ~ 0.29 fm. Both sets of data are 
fitted to V(r) = b + c/r + ar, using the points r = 1,2, y/8, which therefore show no 
deviation. For the remaining points, it is clear that mean-link Tl (circles) shows smaller 
rotational symmetry violations than plaquette TL See also table ^. 
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